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NUMERICAL  MODELING  OF  AURORA 
Volume  I of  the  Calendar  Year  1974 
Annual  Report  to  the  Defense  Nuclear  Agency 

HAES  Report  No.  16 
Section  1 
INTRODUCTION 

A defensive  policy,  particularly  one  of  limited  retaliation, 
requires  a thorough  knowledge  and  understanding  of  the  environ- 
ment within  which  this  policy  may  have  to  be  implemented.  An 
a priori  knowledge  of  the  environment  is  required  to  successfully 
define  and  deploy  the  intelligence  gathering  and  communications  links 
upon  wltieh  such  a posture  must  be  based.  As  these  systems  become 
more  complex,  and  as  they  are  required  to  carry  and  impart  more 
detailed  information,  they  become  more  susceptible  to  the  disturbances 
which  would  be  created  in  a nuclear  exchange. 

The  detonation  of  a high-altitude  nuclear  explosion  (HANE), 
and  even  the  detonation  of  lower  altitude  weapons  with  subsequent 
heave  and  lofting  of  weapon  debris,  can  create  large  scale  ionospheric 
disturbances.  Whether  tlie-.e  disturbances  have  as  their  source  prompt 
gamma-ray,  x-ray,  UV,  optical,  or  IR  emission,  or  later  decay 
betas  or  ionized  debris,  the  net  effect  on  upper  atmosphere  constituents 
is  to  ionize  or  excite  these  species.  The  influx  of  energy  can  initiate 
a complex  series  of  chemical  and  physical  reactions  resulting  in 
enhanced  ionization  levels,  increased  UV,  optical  and  IR  radiations, 
and  perturbations  of  existing  magnetic  and  electric  fields.  The  net 
effect  these  perturbations  will  have  on  defensive  systems  must  be 
thoroughly  understood  so  that  our  defensive  posture,  or  our  ability 
to  effectively  retaliate,  will  not  be  seriously  degraded. 


Note.  Manuscript  submitted  October  it,  11)75. 
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Systems  designers  and  analysts  need  to  have  detailed 
knowledge  of  HANE  effects  but  must  face  the  limitations  of  a test 
ban  treaty.  How  can  systems  be  effectively  exercised  or  tested  without 
resorting  io  actual  nuclear  tests0  The  only  approach  available  is  to 
model  or  simulate,  first  the  HANK  peril  rbatiuns,  and  second,  the 
effects  of  these  perturbations  on  proposed  systems.  However,  any 
large  scale  modeling  or  simulation  effort  at  some  point  must  touch 
base  with  realit'1.  Confidence  in  predictions  or  conclusions  based 
on  the  modeling  increases  as  the  validity  of  the  model  is  established 
by  quantitatively  reproducing  known  conditions. 

In  modeling  HANE  phenomenology,  one  logical  place  to 
establish  such  validity  and  confidence,  is  in  the  aurora.  These 
ionospheric  disturbances,  while  not  directly  applicable  to  HANE 
phenomenology,  have  enough  direct  correlation  in  apparent  observables 
to  effectively  help  bolster  confidence  in  our  ability  to  model  HANE 
disturbances.  While  the  energetics  of  the  emissions  or  particles 
creating  the  disturbances  may  be  different  and  the  ionization  and  heating 
levels  considerably  higher  for  a HANE,  modeling  either  phenomenon 
requires  a detailed  understanding  of  (1)  particle  interactions,  (2)  all 
appropriate  chemical  reactions,  and  (31  potential  changes  or  enhance- 
ments in  ionospheric  fields  and  currents. 

These  three  areas  have  been  addressed  by  the  Plasma  Dynamics 
Branch  during  the  past  year.  The  energetics  of  charged  particle  inter- 
actions has  been  approached  from  several  directions.  Electron  energy 
deposition  has  been  treated  in  a simplified  form  by  assuming  a contin- 
uous energy  loss  model  with  no  pitch  angle  diffusion,  and  in  a more 
detailed  transport  equation  approach  which  takes  into  account  pitch 
angle  diffusion  and  allows  for  discrete  large  angle  scattering  events. 

A first  step  in  understanding  the  interactions  of  protons,  and  heavier 
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debris  ions,  with  the  atmosphere  was  taken  with  the  development  of 
a continuous  loss  model  similar  to  that  developed  for  electrons. 

An  auroral  chemistry  code  has  been  developed  to  permit  an 
examination  of  upper  atmospheric  species  ami  reaction  rates.  The 
code  presently  will  accept  an  energy  deposition  profile  from  either 
of  the  above  deposition  models,  and  by  tracking  14  different  species 
is  capable  of  general ing  both  spatial  and  temporal  profiles  of  selected 
species,  ionization  levels,  emission  rates,  etc. 

The  third  area  of  investigation  involved  the  development  of 
auroral  arc  models.  Wo  were  concerned  here  with  developing  an 
understanding  of  the  gross  morphology  of  the  currents,  fields,  and 
plasma  densities  found  in  auroral  arcs  (small  scale  structure  is 
dealt  with  in  Volume  IV  under  plasma  instabilities').  Modeling  the 
weakly  ionized  part  of  a VV  fireball  and  the  late* time  motion  (>5  minutes) 
of  the  nuclear  disturbed  ionosphere  uses  essentially  the  same  relation- 
ships as  those  for  modeling  auroral  arcs.  The  equations  employed 
included  continuity,  momentum,  and  energy  for  electrons,  ions,  and 
neutrals,  and  an  equation  for  the  electrostatic  potential. 

The  work  performed  at  NHL  during  the  past  year  in  all  three 
of  these  areas  is  summarized  in  the  remainder  of  this  volume.  Much 
of  this  work  has  been  previously  reported  in  NHL  Memo  Reports,  at 
several  symposia,  in  technical  journals,  and  at  DNA  sponsored  meet- 
ings. The  principal  contributors  in  each  technical  area  arc  listed 
as  authors  of  the  section  describing  that  work. 


r 


Reproduced  from 
best  available  copy. 


Socl.im  g 


cham. Mii>  PAirnci  r f.nt.iuiy  dfposjtiox 


i’.S 

IldlPIllil- 

1) . . s 

t,.( 

•1  i.m  1 

,)  K. 

h 'f/rrs 

*11, 

■T.  ' ■*. 

Ap-t  • tng 

:hat 

under 

s'  •'•  di 

r..! : ■ ii 

quite  use 

: ..1  h>  *'■.*. 

pr; 

iic  t ill!) 

♦ • 

i.d 

it  * s i : * 

one  ir  u,- ; 

examine 

tue 

s,  an  ee 

it."  l !: 

;a> 

: H P.  i 

stood  :h. 

* .u.  r or  a : 

i: 

"mat- 

• mo 

» C*  t , j j » 

actio;-..-.  >. 

i the  upper  ad 

r.rjsph 

CA  l *J  * 

r 

lei. 

bo’W',a.  . 

: j > 1 1 : ring 

e-:C 

: gelic 

elf  i.i 

■ .O' 

molec.  o 

F ! loo  a 

tmn. 

•.here 

UlU.'t 

i;. 

• u<"ci; 

d; !l!  ii:  d 

. ('■'•■■  s'.an 

dinh 

t.  i iOf 

, u.u  ■• 

V 

i rans 

iei.mis  of  ihc  interactions 


1 o*(  '.ions,  ."tnz.aior  la.  ! • . an  1 Mt:  ,ir  t < ■ am  ,,  «,?  initialed  by 
tftf.e  i.d acu  ••as.  These  same  alto*,  muma  ><  i r ii.  the  nuclear 
environment  (Starfish  electron  patch,  debris  patch,  etc,  ). 

In  this  section  we  address  the  nroblom  of  modeling  the 
deposition  of  energy  in  the  atmosphere  hv  energetic  elec  trons  and 
protons.  These  particles  a e dirc-.-ted  i'  i the  atmosphere  alone 
magnetic  field  lints  and  dope  dl  tip-  i t-ia  rgv  in  collision. il  processes 
with  the  atomic  and  molecular  ,xp.  . . , : e.j(  ni,  Sii.<  e both  the  total 

number  and  relative  densities  d h<;se  .species  are  altitude  dependent, 
the  energy  deposition  profile  is  a!s<  a): mid"  dependent.  The  collisiunal 
processes  are  interactions  which  can  ionize  the  atmosplieric  constituents 
and  impart  kinetic  energies.  They  therefore  can  initiate  a complex 
set  of  chemical  and  transport  pr<  cc-xsc 

Electron  energy  deposition  has  been  approached  from  t'.vo 
directions.  In  the  first,  a simplified  m<  del  is  generated  in  which  the 
discrete  nature  of  electron  inlerac-ions.  is  ignored  and  the  energy 
deposition  is  treated  in  a continuous  maimer.  In  the  second,  a detailed 

t 


t 

i * 

% 

accounting  of  discrete  collisional  processes  is  maintained  in  an  attempt 
to  create  an  "exact”  model  for  electron  energy  deposition.  Both  models 
are  discussed  in  this  section  and  consideration  is  given  to  the  constraints, 

limitation,  and  applicability  of  each  model.  The  section  terminates  with 

i 

a presentation  of  the  model  used  to  calculate  proton  energy  deposition. 

2.1  ELECTRON  DEPOSITION:  Continuous  Slowing  Down  (CSD) 

Approximation 

l 

To  study  the  qualitative,  and  to  approximate  the  quantitative, 
features  of  atmospheric  electron  deposition,  it  is  desirable  to  have  a 
simple,  fast  running  electron  deposition  code.  To  develop  such  a code 
we  have  used  a generalized  continuous  slowing  down  scheme  for  the 
altitude  dependent  energy  deposition  (see  Julienne  [1974]  for  a complete 
treatment).  Since  the  effect  of  pitch  angle  diffusion  is  noi  treated, 
the  predicted  electron  penetration  is  overestimated  and  energy 
deposition  at  higher  altitudes  is  underestimated.  Given  arbitrary  input 
fluxes,  the  code  is  useful  for  comparing  auroral  energetics,  relative 
excitation  rates,  integrated  emission  rates,  and  the  secondary  electron 
fluxes.  It  is  fast  running,  simple,  and  requires  only  a knowledge  of 
the  electron  inelastic  collisional  cross  sections.  The  code  calculates 
the  volume  production  rates  at  each  altitude  for  all  included  processes. 
One  advantage  of  this  code  over  previous  codes  is  that  it  does  not  rely 
on  analytical  fits  to  the  cross  sections.  The  cross  sections  are 
represented  in  numerical  form  by  relatively  few  points  in  a cubic 
spline  interpolation  scheme.  Thus,  experimental  cross  sections  can 
be  readily  used  where  they  are  available. 

The  flux  of  energetic  primary  electrons,  <*>,  electrons 
cm  sec  * ev  * ster  \ incident  on  the  atmosphere  in  an  aurora 
typically  rangt:;  m energy  from  1 to  30  keV.  The  flux  is  assumed 
to  be  isotropic  over  the  downward  hemisphere  (no  magnetic  field 
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effects),  and  the  vertically  incident  flux  is  obtained  by  integrating 
the  vertical  flux  component  (</>  cos  m)  over  azimuthial  and  pitch 
angles  (v  ~ pitch  angle).  This  gives  electrons  per  square 
centimeter,  per  second,  and  per  eV  vertically  incident.  To 
calculate  the  ionization  and  excitation  levels  of  the  atmosphere  both 
the  primary  and  secondary  electrons  must  be  considered.  The 
secondary  electrons,  predominately  with  energies  below  several 
hundred  volts,  tend  to  be  fairly  isotropic  and  to  lose  their  energy 
locally.  The  present  code  uses  a continuous  slowing  down  (CSD) 
approximation  for  both  primary  and  secondary  electrons. 

The  energy  lost,  dE,  by  electrons  of  energy,  E,  traveling 
a distance,  dx,  through  a gas  of  uniform  density,  N,  is 

dE  = N L(E)  dx  . (2.1) 


The  loss  function,  L(E),  is  determined  by  the  cross  section  for  all 
inelastic  collisions 


E-I 


L(E)  = I En'ew 


cr.  (E) 
i a 1 a' 


?/ 

• A 


12 

WS.  (E,  W)  dW 

JOE 


(2.2) 


The  summation  is  over  all  species,  a,  and  over  the  excitation  and 

ionization  states  for  each  species.  In  the  first  sum,  the  W.fl 

represents  the  energies  of  discrete  excitation  including  the  ionization 

threshold.  The  second  sum  represents  the  creation  of  secondary 

electrons  of  energy  W from  the  various  ionization  states.  The  energy 

differential  cross  section  S.  (E.W)  is  normalized  so  that 

ja  ’ 

E-I. 

]a 

o-|q(E)  = / 2 S (E,W)  dW  . 

la  J0  ia 

<t  represents  the  total  ionization  cross  section  into  state  i. 

jar 
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Once  the  electron  flux  is  known,  the  production  rate  of  some 
state  i in  species  a is  simply 


Pj  (Z) 

lor 


- Nrt  (Z)  j -ME)  <r.a(E)  dE 


(2.3) 


W. 


10t 


If  rr.  represents  an  emission  cross  section  for  some  atomic  or  molecular 


line  then  F is  the  emission  rate  of  thr  *ine. 


The  code  uses  a very  simple  scheme  to  deposit  the  primary 
energy  at  any  altiti  le  z.  The  vertically  incident  flux,  <■'< , is  divided 
into  channels  of  energy  En  to  Enf^,  n ~ 1,  2,N  for  suitably  chosen  N. 
The  flux  in  each  channel  h defined  as 


E . + E 
n+ 1 n 


(/>  = 
li 


/ 2 

"'n 


-2  -1 

<ME)  dE  electrons  cm  sec 


(2.4) 


E . + E 
n-1  n 


Tiie  incident  flux  is  then  treated  as  a series  of  n delta 
functions  of  strength  . Each  delta  function  is  deposited  indepen- 
dently. The  energy  continuously  degrades  from  its  initial  value  E^ 

(at  Z = Z ),  to  a value  of  E (Z)  at  altitude  Z.  Using  t.iese  relationships 
o n 

and  appropriate  cross  section  values  one  may  calculate  the  excitation 
rates  from  degraded  primary  electrons  and  the  total  secondary  electron 
production  rate  spectrum,  which  is  converted  to  a secondar  ■*  flux  using 
the  CSD  approximation.  Similarly,  the  tertiary  flux  can  be  obtained. 
The  procedure  can  be  applied  for  any  order  of  electrons.  However, 
since  it  rapidly  converges  it  is  not  necessary  to  go  beyond  quaternary 

electrons. 

% 

The  cross  sections  currently  used  to  calculate  the  loss 
function  are  shown  in  Table  2. 1.  They  will  be  discussed  in  detail  in 
a separate  report.  The  cross  sections  are  of  two  types:  first,  direct 
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TABLE  2.1 


SUMMARY  OF  02>  AND  O CROSS  SECTIONS 


State 

Reference 

N2  Ionization 

Rapp  and  Englander-Golden  [1965];  Schrarn  et  al. 
[1965]. 

a3E  + 1 

u 

Borst  [1972]. 

B3n 

g 

Stanton  and  St.  John  [1969]  ; Chung  and  Lin  [1972]. 

c3  n 

u 

Finn,  Aarts,  and  Doering  [1972];  Aarts  and 
deHeer  [ 1969] . 

w3a 

u 

Chung  and  Lin,  normalized  to  Chutiian  et  al  [1973]. 

a1  n 

g 

Ajello  [ 1970] (peak);  Green  and  Stolarski  [1972] (high 
energv 

b 1n 

u 

Green  and  Stolarski  [1972], 

b !E+ 

U 

Green  and  Stolarski  [1972]. 

Rydbergs 

Green  rmd  Stolarski  [1972]. 

Vibration  j 

Schu!"  [ 1964],  normalized  to  Spence,  Mauer,  and 
Schulz  [1972]. 

02  Ionization 

Rapp  and  Englander-Golden  [1965];  Schrarn  et  al, 
[1965]. 

a ^A 

g 

Linder  and  Schmidt  [197 1] ; Trajmar  et  al.  [1971  ]. 

' b^+ 

Trajmar  et  al.  [ 1971]. 

V 
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State 


Referer.ee 


6. 1 eV  Green  and  Stolarski,  normalized  to  Trajmar, 
et  al.  [1970]. 

3 v - 

B Green  and  Stolarski,  normalized  to  Trajmar,  et  al. 

9.9  eV  Green  and  Stolarski,  normalized  to  Trajmar,  et  al. 

Rydbergs  Green  and  Stolarski 

0 Ionization  Fite  and  Brackmann  [1959]. 

Henry  et  al.  [1969], 

3S  Henry  et  al.  [1969], 

3s  3S  Zipf  [1974] . 

3s  5S  Zipf  [1974]. 

Sum  (A  1=1,  Green  and  Stolarski 
s=0) 

Sum  (As=l)  Green  and  Stolarski 

Sum  (Al=0)  Green  and  Stolarski 

s=0) 


electron  impact  cross  sections  for  important  states  which  can  be  used 
directly  in  calculation  of  the  loss  function;  second,  cross  sections 
which  are  known  for  specific  processes  i.e.,  optical  excitation 
functions  for  particular  emission  lines,  including  the  effect  of 
cascading.  These  cross  sections  are  not  used  in  the  loss  function 
but  do  permit  the  calculation  of  production  rates  for  specific  transitions. 

Figure  2. 1 shows  a calculated  secondary  electron  flux  at 
altitudes  of  200  and  105  kilometers  for  an  incident  10  keV  electron 
beam  with  unit  flux.  The  clotted  lines  show  the  sum  of  up  and  down 
fluxes  from  Banks  et  al.  [ 1974]  for  the  same  conditions.  The  peak 
deposition  for  Banks  et  al.  is  at  105  kilometers  and  is  at  102  kilometers 
for  the  present  calculations.  The  fluxes  for  the  two  calculations  are 
similar  below  100  keV.  Both  calculations  show  that  the  secondary 
flrjc  is  only  weakly  dependent  on  altitude.  The  continuous  slowing 
down  model  appears  to  give  a reasonable  representation  of  the 
secondary  electron  flux  over  a wide  altitude  range. 

2.2  ELECTRON  DEPOSITION:  Transport  Equation  Approach 

In  modeling  electron  energy  deposition  with  this  formulation,  we 
(see  Strickland  et  al.  f 1974]  for  more  details)  attempt  to  account  for  the 
discrete  nature  of  the  collisional  processes,  in  particular,  large  angle 
scattering  is  allowed  and  the  flux  thus  transported  is  accounted  for. 

The  form  of  the  equation  used  is 

>“jr(rE’E’'')=  -*<te.E>'')  + S(7E‘E>>‘)  (2-5) 

where  '!>  is  the  electron  flux  and  S is  the  source  function  with  both 
energy  and  angular  dependence.  The  source  function  is  defined  as 

S(rE,E,M)  =So(rE.Ef*0  + f G(E,E>,  M')<MrE,E,M)dE’  d/P  . (2.6) 
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FLUX  {cm-2  sec-1  ev"1)  PER  UNIT  INCIDENT 


E(eV) 


Fig.  2.1  — I-ow  energy  electron  flux  at  200  and  105  km  from  a unit 
incident  flux  (1  electron  cm-2  sec-1 ) of  10  dev  electrons.  The  solid 
lines  are  from  the  present  calculation  and  the  dotted  lines  are  from 
Banks  et.  al. 
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The  term  Sq  is  the  initial  source  function,  or  driving  term,  in  the 
equation,  G is  the  redistribution  function  accounting  for  ela.-  tic  and 
inelastic  scattering.  The  terms  in  G arc  given  in  Strickland  and 
Kepple  [1974].  The  independent  variables  are  E the  kinetic  energy  of 
the  electron,  n the  cosine  of  the  pitch  angle,  and  r the  optical  depth 
for  an  electron  of  energy  E.  The  optical  depth  in  differential  form 
is  given  by 


n(z) 


.(E) 


E 

k 


VE) 


dz 


(2.7) 


where  o-  is  the  elastic  cross  section,  and  or,  is  the  kth  inelastic 
e k 

cross  section.  The  sum  includes  ionization  cross  sections.  The 
present  formulation  of  optical  depth  assumes  a one-constituent 
atmosphere. 

The  source  function  represents  the  incident  flux  at  the 
upper  boundary  of  the  atmosphere.  Production  of  electrons  within 
the  atmosphere  due  to  this  energy  source  is  accounted  for  in  the 
integral  term  rather  than  in  the  source  term.  Equation  (2.  5)  is 
solved  in  a matrix  formulation.  The  matrix  elements  are  calculated 
by  approximating  4>  quadratic  in  the  logarithm  of  E,  and  linear  in 
/x within  a given  energy/pitch  angle  cell.  An  iterative  procedure  is 
required  to  calculate  backscattered  fluxes.  The  solution  converges 
in  as  few  as  three  or  four  iterations  using  a scheme  which  projects 
the  iterated  fluxes  to  their  final  values. 

The  term  in  the  redistribution  function,  G (Equation  (2.6)), 
accounting  for  elastic  scattering  is  expressed  as 

Ge  = -j-jgy  .e(E’)ME  -E')pe(E,^,/)  , (2.8) 
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where  p is  the  elastic  angular  redistribution  function  normalized 
e 

in  M, 


pe  afE)  J 
e o 


f 


2 7T 


(r 


d u 


(o)  d</>' 


(2.9) 


d it  /dll  is  the  differential  elastic  cross  section  which  is  assumed  to 
e 

take  the  form  of  the  non-relativistic  screened  Rutherford  cross 
sec  tion, 


do- 

e 

d n 


2 4 


z e 


1 

(1  - cos  0 + 


(2.10) 


8 is  the  scattering  angle  and  tj  is  a screening  parameter.  The 
energy  dependence  of  t?  is  taken  from  Berger  et  al.  [ 1970 J . The 
program  currently  allows  for  a cutoff  o beyond  which  scattering 
is  not  permitted.  This  allows  us  to  determine  the  importance  of 
large  angle  scattering  relative  to  multiple  small  angle  scattering 
and  to  make  a more  comprehensive  comparison  with  Fokker-Planck 
results. 


The  Fokker-Planck  equation  can  be  obtained  by  expanding  the 
integrand  in  Equation  (2.6)  in  a Taylor  series.  Before  doing  this,  we 
rewrite  Equation  (2.  5)  in  a form  more  convenient  for  this  expansion: 


(z,  E,/x)  = -n(z)  o-t(E)<Mz,  E,m  ) + n(z)  J a (AE,  E +AE,  8 ) 


•<£{z,E+AE,  jx  +A/u ) d AEd 


(2.11) 


where  z is  the  distance  along  the  field  lines. 

AE  = E'-E, 

A m = d - n , 

2 -1  -1 

and  t t is  the  sum  of  all  differential  cross  sections  in  cm  ev  ster  . 
The  source  term  SQ  has  been  dropped  from  Equation  (2.  5)  foi  the 
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purposes  of  this  discussion.  The  integrand,  <r<t,  is  now  obtained  by 
expanding  about  E and  m to  second  order,  giving: 

= <‘r’0  iE  + 1 A 

O 

1 2 2 2 

+ 2^5(<r,',)iE  + lMViE  (2.12) 

where 

o-’ffe’  = tr  (AE,  E + AE,  6 )<t(z,  E + AE,  m + A^)  and 
^ = ct(aE,E,0)<*>(z,E,m)  • 


Inserting  Equation  (2.12)  into  Equation  (2.11),  we  obtain  the  Fokker- 
Planck  equation: 


dz 


-^<t<AM>  <*><AE>  + ^-^  <^u2> 


+ ^-^^<ae2>  + 


du 

5m^E 


2 . 2 

5/u. 


<(AmAE  ^ 


(2.13) 


where  each  of  the  diffusion  coefficients  <X>  is  given  by 

<X>  = n(z)  J <r(AE,  E,  0)Xd  AE  dO  . (2.14) 

Equation  (2. 13)  can  be  found  in  Walt  et  al.  [1967]  , where  they  have 
dropped  the  last  two  terms  as  being  small  second  order  ones,  A 
complete  derivation  of  the  diffusion  coefficients  for  Equation  (2.13) 
can  be  found  in  Strickland  et  al.  [1974], 
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The  status  of  auroral  electron  energy  deposition  codes  was 
recently  reviewed  in  a DNA  sponsored  meeting  (June  14,  1974)  wh^re 
the  code  work  by  various  contractors  was  compared.  Participants  in 
this  meeting  included  Naval  Research  Laboratory,  Lockheed  Missiles 
and  Space  Co.  (LMSC),  Mission  Research  Corporation  (MRC)  and 
Visidyne.  The  several  contractors  were  asked  to  supply  calculated 
energy  deposition  rates  for  an  incident  flux  with  an  energy  dependence 


of  the  form 


* (E)  = exp  (-E(keV)/lO) 


(2.15) 


This  flux  is  assumed  isotropic  at  a 250  km  upper  boundary  of  a 
CIRA  (1965)  model  atmosphere.  The  calculations  began  at 
50  keV. 

Figure  2.2  shows  several  calculated  energy  deposition  rates 
using  the  above  flux  and  several  assumed  cutoff  angles  in  Equation  (2.10). 
For  these  particular  results,  our  set  of  N„  cross  sections  were  applied 
to  the  sum  of  the  three  atmospheric  constituents.  As  the  cutoff  value 
of  the  scattering  angle  is  decreased,  the  expected  trend  develops, 
i.e.,  deeper  penetration  of  the  electron  flux  and  lower  loss  rates  at 
high  altitudes. 

The  transfer  results  for  0^  = 180°  are  given  in  Figure  2.3 
together  with  three  other  results.  The  LMSC  curve  is  the  Fckker  - 
Planck  result  by  M.  Walt,  the  MRC  curve  is  based  on  a modified 
continuous  slowing  doAvn  approach  (incorporating  pitch  angle  diffusion) 
by  Tarr  et  al.  and  the  BSM  curve  is  a Monte  Carlo  result  by  Berger 
et  al.  [1974].  The  transfer  result  in  this  figure  is  not  identical  to 
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Fig.  2.2  — Calculated  energy  deposition  rates  for  a flux  defined  by  <t>  = exp  (— E (kcv)/10),  with 
selected  limits  cn  allowed  electron  angular  scattering  (equation  2.10).  The  dashed  curves  are 
results  by  Walt,  et.  al  for  the  same  incident  flux 


that  for  0 = 180°  in  the  previous  figure.  The  three  atmospheric 
constituents  were  treated  explicitly  in  this  last  calculation.  There 
is  good  overall  agreement  between  the  various  methods  used  to 
treat  electron  energy  deposition.  Of  the  differences  that  do  exist 
between  our  results  and  the  others,  the  trend  in  Figure  2.2 
suggests  that  they  are  due,  at  least  in  part,  to  our  more  detailed 
treatment  of  angular  scattering. 


ALTITUDE  {km) 


200 


ENERGY  DEPOSITION  RATE  (ev-cm-3-secH) 


Fig.  2.3  — A comparison  of  several  predicted  energy  deposition  rates  for  the  same  in- 
cident flux.  Note  that  there  is  good  agreement  betweei  these  results,  within  a factor 
of  two,  over  all  altitude;  above  90  km. 


2.3 


PROTON  DEPOSITION 


In  addition  to  energetic  electrons,  other  particles,  such 
as  protons,  are  known  to  play  significant  roles  in  the  creation  of 
natural  aurora.  For  example,  hydrogen  emissions  observed  in  auroral 
displays  can  be  explained  by  protons  becoming  excited  hydrogen  atoms 
through  charge  exchange  processes  with  the  atmospheric  molecules. 

The  modeling  of  proton  energy  deposition  can  be  both  simpler  and 
more  complicated  than  that  of  electron  energy  deposition.  As  charged 
particles,  the  protons  and  electrons  are  constrained  to  spiral  around 
magnetic  field  lines  in  the  earth's  atmosphere.  However,  through 
charge  exchange  processes  an  energetic  proton  can  become  an 
excited  hydrogen  atom  and  in  this  uncharged  state  can  readily  cross 
field  lines.  NRL's  goal  in  modeling  proton  energy  deposition  is  to 
develop  a simplified  model  which  will  be  a useful  first  step  in  under- 
standing and  qualitatively  predicting  the  altitude  dependent  volume 
production  rates  of  the  various  ionized  and  excited  states  of  the 
atmospheric  species.  To  i_alculat3  these  rates,  a model  for  the 
energy  loss  and  resulting  slowing  down  of  die  incident  proton  is 
needed.  The  scheme  adopted  by  NRL  is  the  continuous  slowing  down 
approximation  (CSDA)(described  in  detail  in  Rogerson  and  Davis  [1974]). 
Effects  due  to  the  earth's  magnetic  field  are  neglected  and  the  model  is 
one-dimensional. 

Peterson  [1969]  has  shown  that  the  CSDA  works  very  well 
for  electrons  if  their  energy  is  sufficiently  large  that  the  energy 
loss  in  a single  process  is  a small  fraction  of  the  original  energy. 
Another  limitation  required  by  Peterson  is  that  a large  number  of 
processes  with  widely  distributed  thresholds  be  included. 
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In  the  model  generated  here  proton  energies  of  1 keV  or 
larger  are  considered  whereas  the  excitation  and  ionization  processes 
have  thresholds  of  the  order  of  10  eV.  Hence  in  any  single  process 
the  ratio  of  threshold  to  incident  energy  is  small.  In  view  of  this, 
and  because  many  processes  with  a wide  range  of  thresholds  are 
used  in  the  proton  loss  function,  the  CSDA  should  be  a reasonable 
approximation  for  the  altitude  dependent  energy  deposition  of  protons. 
Comments  made  on  the  limited  applicability  of  the  CSDA  for  electrons 
apply  equally  to  the  treatment  here  for  protons,  i.e. , the  code  will 
predict  too  much  penetration  and  will  ignore  backscatter.  However, 
it  is  a relatively  simple  and  fast  running  code  and  will  be  useful  for 
analyzing  and  comparing  the  same  processes  as  the  electron  CSD 
code.  The  effects  of  hydrogen  atoms  created  by  charge  exchange 
are  included  by  using  the  charge  exchange  formalism  due  to 
Allison  [1958],  Secondary  electron  effects  are  treated  by  a scheme 
due  to  Julienne  [1974]  . 

In  the  CSDA  , instead  of  undergoing  discrete  energy  losses, 
the  pioton  losses  are  modeled  in  a continuously  varying  process. 
Energy  loss  is  described  by 


dE  = -N(Z)L(E)dZ  (2.16) 

where  E is  the  proton  energy,  N(Z)  is  the  local  density  of  the  atmosphere 

at  altitude  Z,  and  L (E)  is  the  proton  loss  function  for  energy  loss  per 

P 

unit  of  length  per  molecule  in  the  vertical  direction.  For  simplification 
only  Ng  is  presently  included  as  an  at.  lospherie  constituent.  The  1965 
Jacchia  [1965]  model  atmosphere  is  use  ! in  this  calculation. 

The  basic  proton  (H+)  reactions  with  Ng  that  are  included  in 
the  loss  function  are  direct  ionization  of  N2>  charge  exchange  ioni- 
zation of  N2  (creating  a hydrogen  atom)  excitation  of  N2,  and  elastic 
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scattering.  The  hydrogen  atom  created  by  charge  exchange  processes 
can  undergo  direct  ionization  of  excitation  of  N , electron  stripping 
of  the  hydrogen  atom  in  collision  with  N'  and  elastic  scattering.  It 
is  assumed  that  an  initially  pure  incoming  proton  flux  becomes  a 
charge  equilibrated  mixture  of  II r and  H resulting  from  the  charge 
changing  reactions.  For  a complete  description  of  all  terms  in  the 
proton  loss  function  L^(E)  of  Equation  (2.10)  the  reader  is  referred 
to  Roger  son  and  Davis  [l974j. 

4- 

Figure  2.4  displays  a number  of  total  production  versus 

altitude  curves  calculated  with  the  CSDA.  The  plots  all  have  roughly 

the  same  shape.  The  500  keV  deposition  curve  produces  few  u*  N0 

ions  than  the  10  keV  and  50  keV  curves  at  higher  altitudes  and  crosses 

over  and  exceeds  these  curves  at  lower  altitudes.  This  is  because 

charge  exchange  is  essentially  turned  off  at  500  keV  and  the  total 

+ 

cross  section  for  Ng  production  is  lower  than  at  10  or  50  keV,  As 

4 

expected,  the  altitude  of  peak  creation  decreases  with  increasing 
proton  energy. 

Figure  2.  5 shows  the  contributions  from  the  various  ion 
producing  reactions  to  total  production  for  an  incident  100  keV 
proton  flux.  At  higher  altitudes,  direct  ionization  by  protons  and 
ionization  by  electrons  are  dominant:  at  lower  altitudes,  where  the 
proton  has  lost  energy,  charge  exchange  and  atom  ionization  become 
more  prominent  until  they  actually  dominate  at  the  altitude  of  peak 
Ng+  production. 

Figure  2.6  shows  the  secondarv  electron  flux  at  a number  of 
altitudes  for  the  same  100  keV  proton  flux.  These  curves  have  roughly 
the  same  shape  and  appear  to  be  typical  of  secondary  electron  distri- 
bution (EMG  13).  The  curve  at  105.5  kilometers  represents  the  secondarv 

4- 

electron  flux  generated  near  the  peak  N'  production. 
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200 


VOLUME  PRODUCTION  RATEfCM'JSEC'') 

2.4  — Total  volume  production  rate  of  N2+  for  selected  incident  proton  energies, 
expected,  the  altitude  of  peak  N2  creation  decreases  with  increasing  proton  energy. 


VOLUME  PRODUCTION  RATE  (C M'3  SEC'1 

Fig.  2.5  — Contributions  to  total  N'2+  production  from  various 
ion  producing  reactions  for  an  incident  100  kev  proton  flux 


The  model  and  results  shown  here  are  only  the  first  step  in 
understanding  the  basic  physical  mechanisms  and  energy  partition 
among  the  various  ionization  levels  in  a proton  aurora.  Because  of 
the  simplicity  of  the  model  the  volume  production  rates  are  probably 
underestimated  and  depths  of  penetration  are  probably  overestimated. 
Planned  future  code  improvements  include:  (1)  using  a distribution 
of  proton  energies  rather  than  delta  function  energies;  (2)  adding  02 
and  0 to  the  atmosphere;  and  (3)  introducing  the  earth's  magnetic 
field.  Also  desirable  at  this  time  would  be  a comparison  between 
the  results  of  this  code  and  other  code  work  modeling  proton  energy 
deposition  and  a comparison  between  experimental  data  acquired  on 
events  such  as  PCA’s  witn  the  results  of  this  code.  In  the  near  future, 
the  results  for  proton  deposition  will  be  extrapolated  to  heavier  ions 
such  as  HANE  debris  patch  particles. 


Section  3 

HIE  AURORAL  CHEMISTRY  CODE 

K.  Hyman.  P.S.  Julienne,  and  D.F.  Strobel 

During  the  past  year,  NRL  has  developed  a computer  code  to 
model  the  chemistry  in  the  atmosphere  during  and  subsequent  to  energy 
deposition  by  auroral  electrons  (see  Hyman  [1974]),  The  code  is  cap- 
able of  generating  altitude  profiles  for  14  atmospheric  species  as  a 
function  of  time  and  for  any  given  incident  electron  spectrum.  The 
code  can  predict  volume  and  column  emissions  of  selected  UV,  visible, 
and  IR  spectral  lines. 

AURORA  is  a one-dimensional  code  that  models  the  chemistry 
in  the  upper  atmosphere  after  onset  of  auroral  precipitation.  The 
production  rates  for  auroral  processes  are  calculated  with  the  electron 
deposition  code  described  in  Section  2.1.  The  processes  included  in 
the  code  are  ionization,  dissociation,  and  dissociative  ionization  of 
Ng  and  Og,  and  ionization  of  0,  all  by  electron  impact.  Additionally, 
a number  of  excitation  processes  in  each  of  these  species  are  moni- 
tored. A list  of  the  electron  production  rates  included  in  AURORA  and 
the  references  for  these  rates  are  given  in  Table  3. 1.  Fourteen 
atmospheric  species  are  tracked  by  AURORA.  The  production  rates 
can  be  scaled  up  or  down  to  simulate  aurora  of  varying  strengths  for 
any  given  spectrum  of  incident  electrons,  and  they  can  be  turned  on  or 
off  to  observe  the  transient  response  of  the  atmosphere  to  auroral 
precipitation. 

In  the  altitude  regime  between  100  and  150  kilometers,  where 
most  of  the  energy  deposition  occurs,  7 of  the  14  species  tracked 
have  short  equilibration  time  constants.  Since  we  are  not  primarily 
concerned  with  processes  occurring  on  a short  time  scale  (approx- 
imately 1 second  or  less)  these  species  are  not  integrated  but  are 
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Table  3.1 

ELECTRON  PRODUCTION  RATE 


N, 

Reference 

Vibrational 

Shulz[  1964]  normalized  to  Spence, 
Mauer  and  Shulz  [1972]. 

3 

A u (Vegard-Kaplan) 

Borst  [1972]. 

O 

B 11  (First -positive) 

Stanton,  St.  John  [1969]  joined  to 
Chung  and  Lin[  1972]. 

3 

C nu  (Second-positive) 

Finn,  Aarts,  Doering  [1972]; 
Aarts  and  DeHeer  [1969]. 

W3A 

u 

Shape  of  Chung  and  Lin  [1972]; 
normalized  to  Chutiian,  Cartwright 
and  Trajmar  [ 1973]- 

a1  n (Lyman-Birge-Hopfield) 

Aiello  [1970]  (peak)  and  Green  and 
Stolarsld  [1972]  (high  energy). 

bXn 

U 

Green  and  Stolarski  [ 1972]. 

b'1!  + 
u 

Green  and  Stolarski  [1972]. 

clnu 

Rydbergs 

. 1_  + 

c Xu  J 

Green  and  Stolarski  [1972], 

Green  and  Stolarski  [1972], 

Ionization 

Rapp  and  Englander -Golden  [1965]; 
Schram,  deHeer,  VanderWiel  and 
Kistenaker  [ 1965]- 

Ratio  3914/ionization 

Stanton  and  St.  John  [ 1 9 69  ] ; 
Srivastava  and  Mirza  [ 1 968]- 
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Table  3.1  (Continued) 


0, 

Reference 

a1  A 

g 

Trajmar,  Cartwright  and  Williams 
[ 1971];  Linder  and  Schmidt[  1971]. 

b4  + 

g 

Trajmar,  Cartwright  and  Williams  [1971] 

6. 1 ev  peak 

Green  and  Stolar ski  [1972],  normalized 
to  Trajmar,  Williams  and  Kuppermann 
[1971]. 

3 

B I (Schumann-Runge) 

Green  and  Stolar  ski  [ 1972],  normalized 
to  Trajmar,  Williams  and  Kuppermann 
[1971]. 

9.  9 ev  peak 

Green  and  Stolar  ski  [ 1972],  normalized 
to  Trajmar,  Williams  and  Kuppermann 
[1971]. 

Rydbergs 

Green  and  Stolarski  [ 1972]. 

Ionization 

Rapp  and  Englander -Golden  [ 1965]; 
Schram,  deHeer,  vanderWiel  and 
Kistenaker  [ 1965]. 

Ratio  b^l  /ionization 

Borst  and  Zipf  [1970]. 

0 

Reference 

XD 

Henry,  Burke  and  Sinfailam  [1969]. 

Henry,  Burke  and  Sinfailam  [ 1969], 

3s  3S  (1304) 

| 

Zipf  (private  communication) 

3s5S  (1356) 

Zipf  (private  communication) 

Sum  ( Af  = 1,  As  - 0; 
A?  = 0,  As  = 0) 

Green  and  Stolarski  [1972]. 

Ionization 

i 

Fite  and  Brackmann  [1959]. 
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determined  by  assuming  that  they  are  in  quasi-equilibrium.  The 

3 *f  4.  1 1 

species  with  short  lifetimes  are:  N9(A  X),  N ' , N , 0(  D),  0(  S), 

+ 2 l “ ^ 

0 ( D}  and  CL(b  1).  The  seven  species  which  have  long  lifetimes 

4 2++  1 

and  arc  therefore  integrated  a-e:  N(  S),  N(  D),  O , NO  , 09(a  S) 

+ 

and  02  . 

Table  3.2  lists  the  chemical  rcactic...,  and  rates  used  in 
AURORA.  Also  in  the  code  are  important  reactions  with  thermal 

4*  4*  4- 

electrons:  dissociative  recombination  of  N2  , NO  , and  02  , and 
electron  impact  deexcitation  of  the  metastable  states.  Initial  con- 
ditions for  the  constituents  N2,  02,  and  O can  be  taken  from  any 
of  the  desired  standard  models  such  as  Jacchia  65  or  71  or  Cira  65 
tables.  Initial  conditions  for  the  minor  species  are  based  on 
Strobel's  [1971]  calculations. 

In  addition  to  predicting  atmospheric  species  densities  as  a 
function  of  altitude  the  code  calculates  both  the  volume  emission  and 
column  emission  of  a variety  of  spectral  lines  in  the  UV,  visible,  and 
IR.  AURORA  presently  monitors  the  3914,  4278,  and  4709  Angstrom 
emissions  from  the  N2+(B£)  state.  These  are  a direct  measure  of  the 
instantaneous  deposition  of  auroral  electrons.  The  McClatchey  [1972] 
atmospheric  absorption  model  for  the  visible  and  infrared  lines  is  used 
to  determine  how  much  of  the  emission  will  survive  the  lower  atmosphere. 

To  exercise  the  code,  and  to  examine  the  buildup  and  decay  of 

various  ionization  levels,  the  electron  spectrum  found  in  Figure  3. 1 

was  used  to  obtain  the  results  presented  below.  This  spectrum  is  a 

slightly  modified  ICECAP  spectrum  (J.  Ulwick,  AFCRL,  private 

communication)  and  corresponds  in  intensity  to  an  ICB  II  aurora. 

The  flux  is  isotropic  downward  (along  the  field)  at  200  kilometers 
• 0 

and  the  magnetic  dip  angle  is  77  . In  the  results  to  follow,  wc  have 
assumed  this  electron  flux  incident  from  time  t = 0 to  t = 10 
seconds.  Major  atmospheric  constituents  were  obtained  from  the 
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•unit'd  input  electron  spectrum  used  to  exercise  AURORA, 
spectrum  is  a slightly  m:  ' fied  ICECAP  spectrum. 


Tabic  3.2 

CHEMISTRY  RATES  IN  AURORA  (Temperature  in  ev) 


Reaction 


(l-C) 

N 

+ O-N  + O 

(2-C) 

N+ 

+ NO— N + NO+ 

(3-C) 

N+ 

+ O0  ~N(2D)  + 0, 

(4-C) 

0+ 

+ 02-0  + 02 

(5-C) 

0+( 

;2d)  + n2-o  + n2 

(6-C) 

0+(2D)  + 02-0  + 02 

(7-C) 

N2+ 

+ °2~  N2  + °2+ 

(8-C) 

N2+ 

+ o2-n2  + o2 

(9-C) 

°2+ 

+ N(2D)-02  + N 

(10-C) 

°2+ 

+ no~o2  + no" 

(1-RN)  N + 

NO-N2  + 0 

(2-RN)  N - 

00—  NO  + 0 

1. 0 x 10 

8.0  x 10' 

2.8  x 10' 

4.8  x 10' 

4.  6 x 10 

1.0  x 10' 

7.0  x 10' 

3.0  x 10' 

3.3  x 10' 

2.9  x 10' 

3.9  x 10' 

2.5x10 

6.3  x 10' 
2.7  x 10' 


Reference 


12  Estimate 

'10  Fehsenfeld  [1972] 

'10  T.  < . 4 McFarland  [1973] 

■10  T (.57)  T ><4 

l l 

'12  T.  ~A  Tj<  . 1 55  McFarland  [1973] 

'10  T.1,2  T.  >.155 

i l 

10  Rutherford  [1971] 

^ Stebbings  [ 1 966] 

'10  Fehsenfeld  [1972] 


12  Ti"'8  Ti  < • 3 McFarland  [1973] 

• 11  14 

11  T.  T.  >.3 

l i 


Strobel [1974] 
Fehsenfeld  [ 1970] 
Garvin  [1973] 


1.3x10  T exp(-.271  Garvin  [1973] 

n j \ 

n' 


t 
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Table  3.2  (Con’t) 


Tabic  3.  2 (Ccm't) 


Reaction 


(8-Q)  N2(A3E^  + N— N2  t- N 5.0x10 


*-*'  2 2 2 

(13-Q)02(a^.i ) 4 0^02  + O 1.0x10 

(14-Q)C2(a1A)  - 02-O2  + 02  4.  1 x 10 

(15-QK)2(b!D  * N2^G2^  n2  2,0x  10 
(le-QlO^lb1^)  • 0„-09  + 09  1.5  x10 


Rate 

Reference 

3,4x.0-9T1/2 

n 

Garvin  [1973] 

4.3  x 10*l2e.xp(-.073  Tn) 

Garvin [1973] 

5.  0 x 10"U 

Young  [ 1968] 

1.  0 x io'10 

Shemanskv  [1971] 

CN 

♦H 

1 

o 

X 

o 

ITD 

2.  G x 10'10 

Sh'  manskv  [1971] 

9.  0 x 10'12 

Shemanskv  [1971] 

1.  0 x 10‘17 

Clark  [1969] 

4.1xlO-17Tn(8) 

Garvin [ 1973] 

-15 

2.  0 x 10 

Zipf  ' ' 9 G 9 ] 

1.5  x 10'lG 

Garvin [1973] 

CIRA  1965  model  5 hour  0.  Temperatures  were  assumed  to  be  the 

ambient  neutral  values  given  ’n  this  CIRA  model.  To  examine  the 

effect  of  elevated  temperatures  we  also  present  the  results  of  a "hot 

ion"  model  in  which,  instead  of  ambient  temperatues,  we  assumed  the 

ion  temperatues  tc  be  1 eV  and  the  electron  temperatures  to  be  0. 1 eV. 

Figures  3.2  and  3.3  show  the  species  concentrations  of  O ( S),  O ( D), 

and  NO+.  The  ratio  of  NO+  to  O4  is  shown  in  Figure  3.4,  0+(2D) 

charge  exchanges  very  rapidly  with  N„  and  09  so  that  it  disappears 

^ ^ + 2 

as  soon  as  the  auroral  flux  is  turned  off.  The  0(D)  profile  will 
be  unchanged  with  time  for  a constant  input  flux  because  it  equilibrates 
almost  instantaneously  and  because  its  production  and  destruction  rates 
are  dependent  on  O and  on  Ng  and  0^  respectively.  These  species 
maintain  their  ambient  values. 

The  production  of  NO+  is  due  primarily  to  the  N^,  O 
rearrangement  reaction  and  O^,  NO  charge  exchange.  NO+  is  destroy- 
ed only  by  dissociative  recombination.  Since  0^  and  disappear 
rapidly  when  the  flux  is  turned  off,  so  does  NO+  production.  Because 
NO+  is  the  terminal  ion,  its  production  rate  is  directly  proportional  to 
the  incident  flux.  In  Figure  3.4,  the  ratio  NO^O^  rises  steadily  after 
the  source  is  turned  off.  This  is  a manifestation  of  the  fact  that  NO+is 
the  terminal  ion  and  disappears  more  slowly,  i.e.,  the  NO+  recombination 
time  exceeds  the  0^  charge  exchange  time  constant.  To  examine  other 
species  and  levels  and  to  consider  the  effect  of  long-term  energy  depo- 

4 

sition  (to  10  seconds),  the  reader  is  referiod  to  Hyman  and  Julienne [ 1974]. 

The  results  presented  here  agree  well  with  the  number 
densities  and  chemical  species  typically  measured  in  IBC  II  auroras. 

It  is  necessary  now  to  obtain  measured  electron  fluxes  and  species 
emissions  to  compare  experimen  al  values  with  results  predicted 
by  AURORA.  Such  an  effort  is  c irrentiy  underway  and  the  results 
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of  this  study  should  be  reported  early  next  year.  Zipf's  [1969]  obscr- 

10  -3 

vations  of  aurora  with  very  high  NO  densities  (greater  than  10  cm  ) 
and  very  large  NO"VO^"  ratios  have  provoked  extreme  interest  in 
examining  the  efficiency  of  converting  the  energy  of  precipitating 
electrons  to  the  formation  of  the  NO  molecules.  In  preliminary 
calculations  we  have  been  able,  by  making  various  speculative 
assumptions,  such  as  increased  ion  temperatures,  to  duplicate  the 
emissions  observed  by  Zipf.  Future  work  will  continue  these 
calculations  and  will  couple  predictions  from  AURORA  with  the 
auroral  arc  modeling  reported  in  the  next  section  to  examine  the 
effects  of  Joule  heating  and  transport  of  species  resulting  from  electric 
fields  and  currents. 
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Section  4 

AURORA  ARC  MODELING 
J.A.  Fedder  and  J.H.  Orens 


In  developing  an  ionospheric  model  to  predict  the  current, 
electric  fields,  and  plasma  densities  in  auroral  arcs,  we  are  working 
in  an  area  which  has  almost  direct  applicability  to  HANE  phenomen- 
ology. Note  that  we  are  not  concerned  here  with  the  high  beta  flows 
that  are  best  handled  by  the  MHD  HANE  codes  (MICE  and  MRHYDE),  but 
with  the  lower  beta  flows  characteristic  of  the  weakly  ionized  part  of 
the  UV  fireball,  of  the  late-time  motion  of  the  nuclear  disturbed 
ionosphere,  and  the  auroral  plasma.  The  correlations  are 
apparent,  e.  g. , the  space  and  time  scales  are  similar  (tens  to 
thousands  of  kilometers  and  minutes  to  hours),  the  boundary  conditions 

are  essentially  the  same,  the  electron  densities  are  similar  (10^ cm 
7 -3 

to  10  cm  ),  and  the  magnetic  field  is  constant  (low  beta). 


The  auroral  arc  modeling  program  at  NRL  started  with  the 
development  of  a simple  two-dimensional  model.  The  results  of  calcu 
iations  made  with  the  model  were  shown  to  be  in  qualitative  agreement 
with  certain  measurements  of  auroral  currents  and  electric  fields. 
This  model  was  then  dimensionally  extended  and  applied  to  the  exam- 
ination of  the  convection  of  the  ionospheric  plasma  in  the  midnight 
region  of  the  auroral  oval.  Once  again,  the  results  obtained  with  this 
model  were  shown  to  be  in  good  qualitative  agreement  with  patterns 
inferred  from  measurements.  Both  modeling  efforts  are  described 
below. 


4. 1 THE  2-D  AURORAL  ARC  MODEL 

The  simple,  two-dimensional  auroral  arc  model  developed 
the  past  year  (see  Fedder  [1974])  was  initiated  to  calculate  auroral  arc 
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currents,  electric  fields,  and  plasma  densities  in  the  ionosphere.  The 
are  model  is  based  on  two  assumptions: 

e Currents  in  the  ionosphere  and  in  the  auroral  are  are 
determined  by  the  conductivity  and  the  ionospheric 
eleetrie  field.  Owing  to  the  inhomogeneous  plasma 
density  and  electric  field  gradients.  Bi/kelano  currents 
flow  parallel  to  the  geomagnetic  field.  The  closure  of  the 
Birkeland  currents  in  the  magnetosphere  is  not  specified. 

• The  ionospheric  electric  ueid  is  the  superposition  of  the 
magnetospherie  convection  eleetrie  field  and  a magnetic 
field-aligned  potential  created  by  the  existence  of  the 
Birkeland  currents  and  a field -aligned  resistivity. 

The  are  model  is  constructed  in  a Cartesian  system  with  the  z-axis 
positive  upward,  the  x-axis  normal  to  the  auroral  are,  and  the  y-axis 
parallel  to  the  are.  The  geomagnetic  field  is  along  the  z-axis,  and 
the  imposed  uniform  magnetospherie  eleetrie  field  is  either  normal 
or  parallel  to  the  are.  Assuming  that  the  vertical  motion  of  ions 
is  unimportant,  that  the  difference  in  number  density  between  ions 
and  electrons  is  negligible,  and  that  the  auroral  arc  and  ambient 
ionosphere  are  ereated  by  an  assumed  ionization  production  rate, 
P(x,z),  the  are  model  is  defined  by  the  following  set  of  equations: 
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dz 
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where  n 
v 

x 
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d j |(  (X,  Z) 
d Z 


o 

ion  density 

x -component  of  ion  velocity 
recombination  rate 
Pederson  and  Hall  mobilities 


E 

x 
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j 


ll 


R(x) 
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II 


(x) 


ionospheric  electric  fields  in  x and  y directions 
field -aligned  Birkeland  current  carried  by  electrons 
ionospheric  current  in  the  x-direction 
field  aligned  resistivity 
length  cf  resistive  region 

Birkeland  current  density  at  top  of  ionosphere,  and 

lower  and  upper  ionospheric  boundaries. 

geomagnetic  field  strength 
speed  of  light. 


These  equations  form  a complete  set  describing  an  auroral  arc  in  the 

ionosphere.  They  can  be  solved  explicitly  in  time  to  arrive  at  a steady- 

state  solution.  Appropriate  descriptions  of  P(x,  z),  a,  R(x),  L,  and 

E and  E must  be  supplied. 
xo  yo 

There  are  two  questionable  assumptions  inherent  in  the  above 
equations:  first,  that  the  ionospheric  electric  field  does  not  vary  with 
height;  and  second  that  the  ions  do  not  move  vertically  For  this 
preliminary  effort,  these  assumptions  were  justified  solely  on  the 
grounds  that  they  simplified  the  numerical  calculations.  As  this  work 
is  continuing,  these  assumptions  are  being  removed. 
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To  examine  this  model  and  make  qualitative  comparisons  with 
measurements,  we  will  assume  that  the  ionization  production  rate 
consists  of  an  ambient  ionization  rate,  and  an  arc  ionization  rate. 


P(,<’z)  = Parab(z) 


P (x.z) 
arc 


We  assume: 


Pamb(z) 


arc 


-3  -1 

3(h  cm  sec  , and 

7.  5 x 10^  exp  [-0.  004  (x  -x^)^  ] 


2 

/ Az  + Bz  + C;  z <120  km 

X 

' exp  [0.  025  { J20 -z)];  z>120km 

where  x and  z are  in  km,  and  the  arc  is  centered  at  x^.  A,  B,  and  C 
are  chosen  so  that 


2 

Az 

+ Bz  + C - 1 

at  z = 120 

Az2 

+ Bz  + C = 0 

at  z = 100 

2Az 

+ B = 0 

at  z = 115 

The  arc  ionization  profile  approximates  the  results  of  Romick  and  Belon 
[1967,  Figure  13],  and  the  ambient  profile  is  a simple  approximation 
based  only  on  expected  levels  and  trends. 

Specifying  a constant  ionospheric  convection  velocity,  we  can 

examine  and  review  the  applicability  of  the  model  in  more  detail.  Assum 

ing  E = 25  m V m -1,  E = 0 , and  a stationary  auroral  arc  source, 
xo  yo 
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will  generate  the  results  shown  in  Figures  4.  1 and  4.  2.  This  model 
should  be  applicable  to  cast -west  aligned  arcs  in  the  evening  sector 
of  the  auroral  oval  with  the  electric  field  and  x-axis  directed  north- 
ward. It  should  also  be  applicable  to  cast-west  aligned  arcs  in  the 
morning  sector,  but  with  the  electric  field  and  x-axis  directed  south- 
ward. 

The  Birkeland  currents  plotted  in  Figure  4.  1 form  oppositely 
directed  sheet  currents  on  either  side  of  the  arc.  They  close  the  ion 
Pederson  current  which  flows  through  the  arc  in  the  direction  of 
increasing  x.  The  ionospheric  electric  field  is  greater  than  the  mag- 
netospheric  field  on  the  edge  of  the  arc,  but  is  50  percent  less  in  the 
center  of  the  arc.  The  "M"  shaped  field  in  the  vicinity  of  the  arc  is 
characteris  ic  of  all  computations  to  date. 

Figure  4.2  shows  computed  plasma  and  current  density  profiles. 
Note  the  distortion  of  the  plasma  contours  at  altitudes  where  the 
Pederson  currents  occur.  In  the  absence  of  an  electric  field,  these 
contours  would  be  tear-drop  shaped.  The  y-directed  ionospheric 
currents  form  an  arc -associated  electrojet  which  is  broadened  by  the 
enhanced  field  at  the  arc  edges. 

In  a recent  article,  Maynard  et  al. f 1 973]  discusses  the  data  for 
electric  fields  in  precipitation  regions.  Allowing  for  a number  of 
contrary  results,  they  conclude  that  the  electric  field  magnitude  and 
particle  precipitation  are  anticorrelated.  This  is  in  agreement  with 
the  above  results,  and  in  fact  it  is  expt  'd  that  more  measurements 
would  agree  with  these  results.  This  is  because  arcs  arc  generally 
observed  to  be  aligned  in  an  cast-west  direction  around  the  auroral 
oval. 
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- Calculated  auroral  electron  Birkelaml  currents  and 
ionospheric  electric  field  normal  to  an  arc 


The  results  for  currents  computed  with  this  model  also  agree 
qualitatively  with  measurements  of  arc  currents.  In  particular,  the 
example  above  shows  an  eastward  elcctrojet  and  the  Birkeland 
current  out  of  the  ionosphere  to  the  north  of  the  current  in  agreement 
with  Cloutier  [1973].  The  fact  that  the  current  density  distribution 
does  not  agree  with  that  which  has  been  measured  is  not  a serious 
fault  of  the  model.  The  distribution  is  a sensitive  function  of  the 
assumed  form  of  the  magnetospheric  convection  field  and  of  the 
anomalous  resistance  which  are  both  unknown  for  the  current  measure- 
ments. In  future  work  we  will  show  how  changes  in  the  assumed  form 
of  the  field  aligned  resistance  can  result  in  different  Birkeland  current 
distributions. 

4.  2 EXTENDED  AURORAL  ARC  MODELING 

As  a next  step  in  improving  both  modeling  capability  and  our 
understanding  of  the  auroral  process,  a numerical  simulation  of  the 
convection  of  ionospheric  plasma  in  the  midnight  region  of  the  auroral 
oval  was  performed.  To  explain  the  general  features  of  the  polar 
region  ionospheric  currents  which  lead  to  magnetic  disturbances, 
Axford  and  Hines  [1961]  introduced  the  idea  of  a magnetospheric 
convection.  There  have  been  numerous  other  papers  since  then  about 
the  nature  and  causes  of  magnetospheric  convection,  a large  number  of 
which  calculate  plasma  convective  flow  patterns  in  the  magnetosphere 
and  ionosphere  using  relatively  simple  models.  Although  these  models 
show  some  agreement  with  data,  they  are  not  completely  convincing. 
This  is  primarily  because  they  are  cither  essentially  one -dimensional 
calculations  or  the  assumptions  on  which  they  are  founded  are  question- 
able. Our  extension  of  some  of  these  ideas  into  multiple  dimensions 
has  resulted  in  excellent  agreement  with  data. 
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In  performing  this  simulation,  the  ionospheric  plasma  continuity 
equation,  the  ionospheric  current  continuity  equation,  and  an  equ  lion 
for  the  ionospheric  electric  potential  must  all  be  solved  simultaneously. 
The  calculations  are  performed  in  a three-dimensional  right-hand 
Cartesian  grid  centered  at  midnight  (x  = 0,  y = 0).  The  x-axis  is 
directed  southward,  the  y eastward,  and  the  calculations  cover  an 
area  1500  km  square  in  x and  y and  105  km  high  (95 < x < 200  km). 

The  simulation  requires  assumed  ionospheric  plasma  production  and 
recombination  rates,  the  magnetospheric  convection  electric  potential, 
and  a resistance  to  the  field  aligned  Birkeland  current.  The  mathematical 
formulation  extends  the  previously  described  auroral  arc  equations  into 
a third  dimension. 

The  ion  continuity  equation  in  the  auroral  ionospheric  E -region 
is  def  > ied  as 


dn 

ei 


dx 


n v 

x 


o 2 
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where 

n = ion  density 

P = ion  production  rate 
a = volume  recombination  rate,  and 

v , v , the  ion  velocities,  are  defined  as 

x y 
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Vy  " MP  + Mh  "B“ 
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where 


- Pedersen  and  Hall  mobilities 

- speed  of  light 

= geomagnetic  field  strength  (=0.  5 Gauss),  and 

= x and  y components  of  ionospheric  convection 
electric  field. 


Diffusive  transport  and  vertical  transport  of  ionization  have  been 
neglected  since  it  is  expected  that  horizontal  convective  transport 
dominates. 


The  ionospheric  current  continuity  equation  is 


where 


where 


j 

d z 


5 j 


<?x  Jx  9 y ^y 


j = Birkeland  current  density  (assumed  carried 
by  electrons) 

j^,  j = horizontal  ionospheric  currents 
qn  (vx  " vx  qn  (vy  " VyG) 


q = electron  unit  charge,  and 
6 6 

v v = electron  velocities 
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For  the  electric  field, 


E - - V 


* = * + J..  R 

m H 
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where 

= ionosphere  electric  potential, 

^ -assumed  magnetosphere  potential, 

J (!  total  Birkeland  current  entering  and  leaving  iono- 
sphere, and 

R - assumed  anomalous  resistance  to  J(| 

Figure  4.  3 shows  the  assumed  magnetos])hcric  convection  pattern 
as  it  maps  along  field  lines  onto  the  ionosphere.  The  resulting 
electric  field  in  the  polar  cap  ionosphere  is  25  niv  m westward  at 
midnight. 

The  ionospheric  convection  pattern  is  shown  in  Figure  4.4, 

It  is  considerably  different  from  the  magnetospheric  pattern  shewn 
in  Figure  4.  3 . The  difference  is  caused  by  a field  aligned  potential 
between  the  magnetosphere  and  ionosphere.  The  most  important 
feature  of  the  pattern  is  the  westward  shift  of  the  east -west  convective 
discontinuity  before  midnight.  This  behavior  has  been  reported  by 
a number  of  experimenters.  In  particular,  one  should  note  the 
similarity  between  the  calculated  pattern  and  the  pattern  inferred 
from  measurements  by  Heppner  j 1 972 ]. 

The  southward  electric  field  along  three  north-south  planes 
through  the  auroral  oval  is  shown  in  Figure  4.  5.  The  changing  field 
direction  is  very  similar  to  a measurement  of  Gurnett  and  Frank  [1973], 
The  Birkeland  current  density  patterns  a*  ' presented  in  Figure  4.  G. 

The  pattern  seen  here  is  the  same  as  that  observed  in  the  measurements 
of  Armstrong  and  Zmuda  [1974],  The  very  large  amplitude  of  the 
current  out  of  the  ionosphere  at  the  northern  border  west  of  midnight 
is  largely  responsible  for  closing  the  eastward  and  westward  electrojets. 


IONOSPHERE  ELECTRIC  POTENTIAL 


Kin.  4.  t — The  ionospheric  convection  pattern  which  results  from  polarization  of  the 
auroral  oval  and  the  field-aliened  potential  drops.  Note  the  multiple  reversals  to  the 
west  of  midnight. 


Figure  4.7  is  a contour  map  of  the  auroral  electrojet  current. 

The  most  important  feature  here  is  the  extension,  some  500  km  to  the 
west  of  midnight,  of  the  westward  electrojet.  This  is  a feature  of  the 
westward  jet  that  has  been,  recognized  from  geomagnetic  observations 
by  Harang  [1945]  and  Heppner  [1972].  The  overlap  of  the  eastward  and 
westward  electrojets  in  the  southern  part  of  the  oval  has  also  been 
observed  in  geomagnetic  data. 

Based  upon  the  results  generated  with  the  above  model,  we  may 
make  some  tentative  conclusions  concerning  certain  features  of  the 
auroral  oval.  The  Birkeland  current  systems  appear  to  close  a sizable 
portion  of  the  intense  currents  in  the  midnight  oval.  This  includes  not 
only  the  north-south  Pedersen  currents,  but  also  the  eastward  and 
westward  electrojets.  The  Birkeland  currents  occur  in  sheet-like 
structures  at  the  northern  and  southern  edges  of  the  oval  for  two 
reasons;  first  because  of  the  imposed  magnetospheric  electric  field 
pattern,  and  second  because  of  the  essentially  east-west  alignment 
of  the  auroral  plasma  density  enhancement. 

With  satellite  photography  of  the  aurora  [Snyder  et  al.,  1974],  it 
has  become  clear  that  the  aurora  occurs  in  two  forms;  diffuse 
aurora  and  discrete  aurora  [Akasofu.  1974].  Comparing  the 
patterns  of  these  two  forms  with  the  Birkeland  current  patterns 
shown  abo'v'e,  indicates  that  the  discrete  aurora  occurs  primarily  in 
the  region  of  large  amplitude  upward  Birkeland  currents  and  field 
aligned  potential,  while  diffuse  aurora  appear  to  occur  in  regions 
where  the  field -aligned  potential  either  opposes  electron  precipitation 
or  is  weak. 


Section  f> 


SUMMARY  AND  CONCLUSIONS 

In  the  introduction,  we  briefly  discuss e.  'he  rationale  for 
initiating  efforts  to  numerically  simulate  auroral  processes.  The 
impulsive,  transient  disturbances  of  the  ionosphere  anti  atmosphere 
found  in  an  aurora  are  not  as  great  as  those  found  near  a HAKE, 
but  are  similar  to  those  reached  in  the  wings  oi  a I’ V fireball, 
in  a beta  batch,  or  ai  late  tunes  in  atmospheric  heave.  Many  of 
the  mathematical  relationships,  necessary  to  the  simulation  of  BANE 
phenomenology  are  also  used  in  simulating  auroral  processes.  More 
importantly,  many  of  these  relationsmps  are  not  amenable  tc  direct 
analytic,  solution  and  must  therefore  be  solved  numerically. 

Since  we  are  unable  to  conduct  BANE  tests,  and  because 
much  of  the  data  acquired  during  previous  HAKE  tests  are  not 
applicable  to  the  technical  questions  raised  today,  we  must  rely 
heavily  on  our  ability  to  simulate  HAKE  phenomenology'  in  the  paper 
exercise  of  certain  defensive  systems  and  strategies.  Confidence 
in  the  results  of  these  exercises  can  best  be  established  when  as 
much  of  the  overall  simulation  model  as  possible  has  been  exercised 
to  reproduce  appropriate  natural  or  man-made  atmospheric  distur- 
bances. The  largest  available  natural  disturbance  against  which  the 
simulation  model  can  be  checked,  one  which  will  regularly  occur  for 
repeated  checks  and  investigation,  is  the  aurora. 

In  developing  an  auroral  model,  there  are  three  technical 
areas  which  must  be  addressed:  first  is  the  deposition  of  energy 
into  the  upper  atmosphere  bv  energetic  electrons  and  protons,  second 
is -the  complex  series  of  chemical  reactions  initiated  by  this  energy’ 
deposition,  and  third  is  the  perturbation  or  creation  of  fields  and 
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currents  in  the  ionosphere.  All  three  areas  have  direct  applicability 
to  HANE  modeling;  however,  proper  scaling  relationships  or  modifi- 
cations must  be  used.  Each  of  these  areas  has  been  addressed  by 
the  Plasma  Dynamics  Branch  during  the  past  year  and  significant 
progress  has  been  made  in  understanding  both  auroral  and  HANE 
phenomenology. 

During  the  past  year  we  have  developed  simplified  fast 
running  codes  to  model  electron  and  proton  energy  deposition  in  the 
atmosphere.  Both  codes  use  a continuous  slowing  down  approximation 
for  the  energy'  loss  and  handle  cross-field  energy'  diffusion  in  a mini- 
mal fashion.  They  therefore  tend  to  underestimate  energy’  deposition 
at  higher  altitudes  and  to  overestimate  atmospheric  penetration.  The 
codes  are  extremely  useful  however  for  making  qualitative  comparisons 
between  calculated  energy  deposition,  ionization  rates,  and  optical 
emissions  for  different  input  spectra  and  conditions. 

Electron  energy  deposition  has  also  been  modeled  in  more 
detail  using  a transport  formulation  which  attempts  to  account  for 
discrete  energy  losses  and  for  large  angle  scattering,  and  which 
allows  many  angles  to  be  included  in  specifying  the  input  flux..  A 
one -constituent  Ng  atmosphere  is  currently  assumed,  but  0 and 
are  now  being  added. 

Early  results  from  his  detailed  model  were  presented  at 
a June  1974  DNA  meeting.  At  that  time  results  from  this  deposition 
code  were  shown  to  predict  higher  backscatter,  shallower  penetration, 
and  higher  deposition  at  high  altitudes,  than  either  the  MRC  or  LMSC 
codes  (its  behavior  relative  to  the  CSD  approach  is  as  expected--see 
Figure  2.2).  The  model  was  further  validated  by  confining  scatter- 
ing events  to  small  angles  to  approximate  the  Fokker-Planck  solutions 
and  reproduce  the  reported  LMSC  energy  deposition  profile. 
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Energy  depo-ition  studies  \ ill  continue  in  several  areas 
over  the  next  year,  i’i'.e  transport  code  will  be  exercised  to  develop 
a better  understanding  of  elect' on  energy  oss  processes,  and  in 
particular  of  the  magnbuce  ‘if  predicted  backscauer.  The  high 
backscatter  predicted  ny  tins  code  is  a result  of  the  basically 
differ enr  vac.  a '.a,  ti  * energy  ics.->  process  is  modeled  and  it 
is  therefore  ju-cesr  avy  to  obtain  a del  a*  led  understanding  of  the 
role  which  this  .ntoriVTion  plays  in  : ’ne  energy  loss  process.  This 
same  rode  viU  be  tppl'ed  fo  tlm  Starfish  detonation  in  an  attempt 
to  model  the  observed  emissions  from  the  ”hct"  electron  patch. 

The  proton  eneigy  deposition  work  will  be  continued  to  model  selected 
PCA  event  data,  to  handle  the  energy  loss  process  in  a more  realistic 
manner  and  to  extend  the  model  to  heavier  particles  such  as  debris 
ions. 


The  auroral  chemistry  code  developed  during  the  past  year 
will  accept  energy  deposition  rates  from  the  above  models  and  follow 
the  many  reaction  processes  initiated  by  such  a deposition  to  predict 
the  ionization  levels,  and  emission  rates  initiated.  The  deposition 
rates  from  the  above  codes  can  be  scaled  to  simulate  aurora  of  vary- 
ing strengths  for  any  selected  input  spectrum.  The  transient  response 
of  the  atmosphere  may  be  observed  by  turning  the  deposition  rates  on 
and  off. 


The  processes  incorporated  have  been  described  previously, 
but  the  more  important  ones  included  were: 

2 

• Dissociation  of  N2  as  a source  of  NCD)  which,  reacting 
with  Og,  is  the  dominant  source  of  NO. 

• The  decay  of  certain  highly  excited  states  of  Ng  to  supply 
a strong  UV  source.  This  UV  ionizes  Og  and  O and 
dissociates  O2. 
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• Thermal  electron  reactions  sucn  as  dissociative 

recombination  of  N*.  NO+,  and  and  electron  impact 
deexcitation  of  metastable  states. 

Using  a scaled  ICECAP  spectrum  obtained  from  AFCRL  measurements, 
the  code  has  been  exercised  to  demonstrate  the  predictions  obtainable. 
The  qualitative  features  of  the  predictions  were  as  expected.  They  are 
now  being  validated  quantitatively. 

One  phase  of  this  validation  will  test  both  the  electron 
deposition  code  and  the  auroral  chemistry  code.  Measurements  were 
made  during  ICECAP  72  of  primary  electron  flux,  electron-ion 
densities,  positive  ion  densities,  positive  ion  concentrations,  and 
various  emission  rates  for  an  1BC  11  aurora.  The  data  was  acquired 
with  an  instrumented  Black  Br°ut  rocket,  on  both  ascent  and  descent, 
that  apparently  passed  just  under  the  lower  edge  of  the  observed  arc. 

By  using  measured  electron  fluxes  in  the  deposition  code  and  then 
generating  emission  rates  with  the  auroral  chemistry  code,  we 
should  be  able  to  reproduce  the  ICECAP  field  measurements. 

Modeling  the  gross  morphology  of  auroral  arcs,  and  extending 
this  work  to  the  midnight  auroral  oval,  has  given  us  confidence  that  we 
can  calculate  large  scale  polarization  and  plasma  density  in  arcs  and 
in  large  ionospheric  plasmas.  The  current  systems  generated  in 
these  models  can  be  used  to  study  the  stability  of  auroral  arcs  to 
perturbations  such  as  gradient  drift  waves.  This  capability  is  also 
necessary  to  provide  background  conditions  for  the  striation  and 
scintillation  work  described  in  Volume  IV. 

Currently,  both  the  auroral  arc  model  and  the  extension  to 
the  midnight  oval  are  steady  state  calculations  (in  the  sense  that 
the  driving  functions  are  constant).  During  the  next  year  this 
modeling  will  be  extended  to  include  time  dependence.  With  this 
improvement  in  the  model,  transient  fields  and  currents  can  be 
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followed,  :uid  the  auroral  response  to  impulsive  inputs  can  be 
examined.  Also  during  the  next  year,  the  arc  model  will  be  coupled 
with  the  auroral  chemistry  code  to  predict  emission  rates  and  ioni- 
zation levels  which  are  both  influent  ed  by  auroral  currents. 

The  since  kuions  performed  to  date  with  the  auroral  model- 
ing codes  all  have  been  qualitatively  verified  against  a number  of 
observations  and  measurements.  This  bolsters  our  confidence  in 
these  models  and  in  our  ability  to  extend  the  techniques  to  HANE 
phenomenology.  A more  detailed  quantitative  check  of  these  same 
codes  is  necessary  as  is  more  predictive  work  on  expected  emission 
rates,  communications  disruptions,  etc. 

The  ongoing  HAES  program,  with  experimental  efforts  such 
as  ICECAP,  EXCEDE,  SPIRE,  etc.,  should  provide  an  excellent 
base  against  which  the  auroral  code  work  developed  during  the  past 
year  may  be  checked.  Also,  other  data  sources,  such  as  measure- 
ments on  PCA  events,  should  be  regularly  reviewed  for  use  in  this 
code  validation.  To  optimize  the  use  of  such  data  it  is  suggested  that, 
where  feasible,  an  iterative  feedback  loop  be  established  between  the 
measurement  program  and  the  modeling  or  simulation  effort.  With 
such  a communications  Link,  the  simulation  effort  can  be  provided  with 
advance  information  on  sensors  to  be  used  on  experimental  plans,  and 
on  operational  parameters  such  as  ejected  EXCEDE  beam  energy  and 
spectrum.  Also,  the  experimentalist  may  receive  valuable  input  from 
the  modeler  on  those  elements  of  his  measurement  program  which  could 
be  improved  to  provide  more  useful  and  usable  data.  Such  a link  could 
also  provide  for  a more  timely  exchange  of  information  so  that  the 
modeler  is  not  continually  working  on  old  data. 
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